Positional, Reorientational and Bond Orientational Order in DNA Mesophases 
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We investigate the orientational order of transverse polarization vectors of long, stiff polymer 
molecules and their coupling to bond orientational and positional order in high density mesophases. 
Homogeneous ordering of transverse polarization vector promotes distortions in the hexatic phase, 
whereas inhomogeneous ordering precipitates crystalization of the 2D sections with different orien- 
tations of the transverse polarization vector on each molecule in the unit cell. We propose possible 
scenarios for going from the hexatic phase, through the distorted hexatic phase to the crystalline 
phase with an orthorhombic unit cell observed experimentally for the case of DNA. 



The microscopic structure of high density DNA 
mesophases is not only of utmost importance for poly- 
mer science in general [|l| but has a fundamental bearing 
also on topics as far removed from one another as the 
nature of the phase diagram of magnetic vortex arrays 
in type-II superconductors or the structure of geno- 
somes used for DNA transfection in gene therapy related 
studies 1^]. By virtue of conventional wisdom at high 
densities DNA should form crystals of hexagonal sym- 
metry in the plane perpendicular to the long axes of the 
molecules, with long range positional order, and nematic 
(or better cholesteric) arrays with short range positional 
order at densities intermediate between the crystalline 
and isotropic phases. There are however a few refine- 
ments one should add to this picture. 

First of all it was found experimentally that DNA in 
fact forms several crystalline phases [Q,^. At very high 
densities Li+-DNA makes a crystalline phase with or- 
thorhombic symmetry that implies a distorted hexagonal 
unit cell perpendicular to the long axes. The probable 
reason for this is that there are angular frustrations be- 
tween neighboring DNA molecules that make hexagonal 
local symmetry energetically costly. A similar situation 
can also be encountered in frustrated spin systems, such 
as in antiferromagnets on a triangular lattice or alkyl- 
chain systems Q . In all these cases distorting the hexag- 
onal equilateral into isosceles triangles could lower the 
angular part in the interaction energy. In this way two 
pairs of molecules are closer to each other, maintaining 
optimal angles, while the third pair is further apart and 
can be in a non-optimal configuration. 

In addition just below the crystalline phase DNA forms 
a line hexatic phase with short range positional order 
and long range bond orientational order Q]. The vari- 
ation of the positional correlation length as a function 
of DNA density shows that positional order within this 
phase is more liquid-like (shorter correlation length) the 
more DNA density is increased This trend is sur- 
prising and counterintuitive: one would expect the DNA 
array to exhibit increasingly longer ranged positional or- 
der approaching the crystalline phase where it becomes 



(ideally) infinite. Again the progressive disordering of 
DNA could be due |^ to increasing angular frustration 
of molecules as they try to satisfy both the positional 
and the angular constraints imposed by the interaction 
potential. 

The question of exactly how angular frustrations could 
affect crystalline and bond-orientational order in DNA 
arrays will be addressed in this contribution. Angu- 
lar frustrations of course correspond to angle dependent 
terms in the interaction free energy between DNAs. Re- 
cent theoretical investigations ^ have indeed made it 
clear that at small enough interaxial separations the in- 
teraction free energy depends crucially on the mutual 
orientation of the two interacting molecules. This orien- 
tation can be specified in any 2D plane perpendicular to 
the long axes, by giving e.g. the position of the major 
groove of one molecule with respect to the line joining a 
pair of them and of the major groove of the other one 
with respect to the first one This amounts to effec- 
tively defining a 2D vector associated with each of the 
interacting molecules (we call it transverse polarization 
vector, p), perpendicular to their long axes, and the in- 
teraction between them will depend on their separation 
as well as their orientation described by their transverse 
polarization vectors. Basing our hypothesis loosely on 
Q we assume that at low DNA density the interaction 
does not depend on the mutual orientation of molecules, 
at intermediate densities the interaction is minimal for 
parallel p orientation while it is minimal for some finite 
angle between ps for larger densities. 

We assume that the molecules are stiff enough so that 
ordering of all 2D sections perpendicular to their long 
axes are the same. We consider consecutively two pos- 
sible situations. First at intermediate densities the pre- 
ferred angle between transverse polarization vectors in 
all 2D sections is assumed to be zero. The transverse 
polarization order is thus homogeneous and we describe 
it with a 2D constant vector p, in the standard com- 
plex form p = \p\e^'^ , where (j) is the local angle be- 
tween the transverse polarization vector and some pre- 
ferred axis. In the isotropic 2D liquid, hexatic order is 



1 



associated with two symmetric combinations of compo- 
nents of a 6th rank tensor. In the complex form the 
hexatic order parameter can be written in terms of the 
angle between molecular bonds in the form i/j — |-0|e*®^, 
where Q < 9 < The free energy of hexatic ordering 
now depends only on J- hex = -^Tiea: ( | V' P ) j while the free 
energy of polar correlations can be written analogously 
as J-poi = J-poi{\p\'^) ■ The coupling of these two types 
of order can be introduced through an interaction term 
^int = -^mt ( I I bl^ cos 6((/) — 6)) ■ In this picture the free 
energy of the isotropic liquid can become unstable with 
respect to hexatic as well as polar order parameters. 

The total free energy T — Thex + ^poi + -^mt can now 
be written as 



+ ci|V'|b|^cos6(0-6i) + 



b2\p\' 



(1) 



This free energy has five different stable solutions as ob- 
tained by standard methods. Since our main interest is 
in the transition from the hexatic to the distorted hex- 
atic phases, it is sufficient to present only a typical cross- 
section of the phase diagram (described in the coordi- 
nates oi, a2 and ci) in the region where the isotropic 
liquid is already unstable, oi < 0. In this cross-section 
with coordinates (a2, C2) the phase diagram has a cusp- 
like form. Fig. |l|. 

The relevant stable solutions of Eq. ^ are in order of 
increasing density: 

1. Isotropic 2D liquid: \ip\ = 0; \p\ = 

2. 2D hexatic phase: |?/'| ^ 0; \p\ = 

3. Three different distorted hexatic phases with the prop- 
erty: IV'I ^ 0; \p\ ^ 0; 



(A) One with: sin 6(0 -6*) = and cos6((/)- 61) = 1, 
thus 4> = 9 + ^(2?T. -I- 1). Obviously here the polar 
vectors are directed along the molecular bonds. 

(B) One with: sin6(0 - 9) ^ and cos6(0 - 9) = -1, 
thus <j) = 9 + Here the polar vectors make an 
angle ^ with the bond directions. 

(C) Besides these two there also exists a lower- 
symmetry phase with a general angle between the 
polar vectors and the bond directions depending on 
external fields, such as density, and varying from 0° 
to 30°. 

The point orientational symmetry in 2D planes of these 
phases is Cqv for the hexatic phase, Cs for distorted hex- 
atic phases (A) and (B) and Ci for the distorted hexatic 
phase (C). 

In a typical case the transitions hexatic — > distorted 
hexatic are of the second order. The transition from 
the hexatic phase to the low symmetry distorted hexatic 
phase (C) can be either the succession of two indirect 
2nd order transitions or one direct 1st order transition, 
see Fig. |l|. The distorted hexatic phases A and B are 
isostructural and therefore at least in principle there ex- 
ists a continuous way to go between them without any 



phase transition. This scenario implies that the hexatic 
order parameter jV'l decreases in these distorted hexatic 
phases until it vanishes, goes through zero and reemerges 
with a different orientation of the polar vector with re- 
spect to the local hexagonal axes. More formally this 
scenario implies that (p ~ 9 should go from 0° to 30°, a 
path which can be continuous only if it crosses the point 
lipl = 0, where 9 is not defined. If we assume that the 
hexatic order is strong enough and vanishes at no point in 
the phase diagram, then the only way from 9+^ (2n -|- 1) 
to 6* + ^ distorted hexatic is through two phase transi- 
tions. 
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FIG. 1. The 2D section of the phase diagram in coordi- 
nates 02 and ci for ai < 0, Eq. ^. The bold lines mark 
second order phase transition lines. The order of the trans- 
verse polarization vector can induce a transition into three 
types of distorted hexatic phases, where the polar ordering is 
coupled to the uniaxial distortions of the hexatic "unit cell" 
(there is no real unit cell in a hexatic). The hexatic order is 
indicated with dotted a hexagon. The directions of the trans- 
verse polarization vector in phase (C) are in between (A) and 
(B) as indicated by the arrows on the figure. The resulting 
elementary unit cells are shown by bold lines 

The deformation of hexatic order (not shown in Fig. |l|) 
of course follows the symmetry changes of the transverse 
polarization vector. It has to be uniaxial with one axis 
pointing in the direction of the transverse polarization 
vectors of the molecules. This can be seen very easily 
if we introduce a parameter quantifying the local distor- 
tion of the bond orientational (hexatic) field e = |e| e*^", 
where e is a combination of the plane strain tensor 
components: e ^ e^x - ^yy + '2i£xy; e* = e^x - Syy -"^iexy 

The coupling terms that do exist are between the trans- 
verse polarization vector and the bond order distortion 
parameter and have the form |p|^|e| cos2((/) — a). This 
term induces the molecular bond distortion as a sec- 
ondary (improper) order parameter, when p is different 
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from zero. In principle we should also add the cou- 
pling term between the distortion parameter and the 
hexatic order parameter that should be of the form 
|e|'^|-!/)| cos 3 (a — 9). This term shows that hexatic order 
in itself can not induce a distortion and a coupling of the 
distortion parameter to the polar order is crucial. 

We now consider the case of larger densities where ac- 
cording to theoretical estimates |9| the preferred angle 
between transverse polarization vectors of neighboring 
molecules is non-zero. The transverse polarization or- 
der is thus inhomogeneous and we describe it with p(p), 
where p — {x,y). In this case the distorted hexatic 
phases can give way to lower symmetry phases charac- 
terized by a condensation of waves of the transverse po- 
larization vector. The condensation of these modes in- 
variably induces also long range positional order pushing 
the system from a distorted hexatic phase into a 2D crys- 
talline phase with an orthorhombic unit cell. 

If the wave condensation starts from the distorted hex- 
atic phase A ( or B ) , then, in general, for transverse po- 
larization modes, there exist two qualitatively different 
choices of the wave vector k directions with respect to 
the symmetry plane of the distorted hexatic order Fig. 
|l|. The direction of the mode can either find itself in the 
symmetry plane, in which case it can be described with 
a single wave vector, or it can point in a general direc- 
tion, when it has to be described with the complete set 
of (two) directional wave vectors. 

The first case can evidently describe only one dimen- 
sional modes. Indeed, here the direction of the mode 
along y-axis, characterised by ky, is completely decou- 
pled from the component along the x-axis. We thus 
obtain order parameter profiles which describe a one- 
dimensional modulation of the transverse polarization 
orientation. Note that the polarization p is already or- 
dered in the homogeneously distorted hexatic phases ; 
then the wave-like reorientation of p in the (x,y)-plain 
is described by the wave of the pseudo-vector perpen- 
dicular to this plain and parallel to long molecular axes. 
The angle r2(p) of the reorientation of p with respect 
to its direction in the homogeneously distorted hexatic 
phase (i.e. A-phase) is given by the one-dimensional wave 
r2(p) — \n\ cos(fcj,y -I- /3), where |r2| is the module and /3 
is the phase of the corresponding order parameter (OP). 

One of the possible realizations of the resulting ordered 
structure is presented as 1 in Fig. 2. Of course, the struc- 
ture should depend on the relation between the wave- 
length of the modulation and the average distance be- 
tween the molecules. This relation cannot be obtained in 
the framework of the phenomenological theory adopted in 
the present work and depends on the constants of the mi- 
croscopic interaction between DNA molecules. For defi- 
niteness, we present here (structure 1 in Fig. 2) the limit 
case, where the wave vector of the modulation is taken to 
be ky = ^Gy. Here b is the long side of the orthorhombic 
"local unit cell" of the homogeneously distorted hexatic 
phase A (Fig. 1) ; is a unit vector along the y-axis. 
Then the angles of the polarization p directions in adja- 



cent layers, taken with respect to the y-axis are |J7|and 
— |r2| respectively. 

The only scalar invariant possible in the free energy 
descibing this transition is |r2p with the free energy itself 
assuming the form 

T = di\n\^ + d2\n\^ + . . . . (2) 

Typically this free energy describes one 2nd order phase 
transition at di = 0. In the x direction, being the di- 
rection perpendicular to the direction of the condensed 
transverse polarization mode, the system remains liquid 
after the phase transition, with continuous translational 
symmetry. 

In the second case the wave vector k is out of the 
symmetry plane of the 2D section of the system. The 
OP describing the transverse polarization reorientation 
mode is in general four dimensional. Mathematically the 
irreducible representation has two components but phys- 
ically it has to have four components (if the waves are 
to be real). Therefore r2i(p) = \ e'^^e'^'P ; r22(p) = 
in^le'P-e^-P ; U^ip) = n\{p) ; and n^{p) = ^^{p) 
, where ki and k2 have the components ki= (k^. ky) 
and \<i2— {—kx.ky) along the symmetry axes of the "unit 
cell" of the homogeneously distorted hexatic phase ; an- 
gles jJliI and |r22l are two, in general independent wave 
amplitudes. The most general form of the resulting 
wave, which again has to be real, is r2(p) — 0,i — 
(|r2i| cos(kip -I- Pi) + |rj2| cos(k2p -I- 132)) , with Pi and 
(52 being the initial phases of the OP components that 
can be made zero by an appropriate choice of the origin. 

The free energy of the transverse polarization reorien- 
tation described by this OP depends only on two indepen- 
dent invariants : Ii = |r2i |^ -I- |J72 1^ and I2 — \^i^ 1^2 
The corresponding free energy thus assumes the form 

T = dih+d2ll + hh + --- (3) 

Apart from the distorted hexatic phases A, B or C 
with homogeneous transverse polarization order jJlij = 
0; IO2I — already dealt with above, the minimization 
of this free energy gives three different stable states: 

2a jJliI ^ 0; |r22| = or another domain where the 
roles of the fields |r2i| and |r22| are reversed, 

2b |r2i| ~ |J^2|, with either Oi = r22 or fZi = —$12, 

2c or the most general form {) ^ 7^ \^2\ 7^0. 

In what follows we shall limit ourselves to the distorted 
hexatic A and the phases with nonhomogcneous trans- 
verse polarization order that are derived from it, since the 
positional order in the crystalline phase of DNA strongly 
favors this type of distortion S]. 

The structure 2a, see Fig. ^has discrete translational 
order only in one dimension with (as opposed to structure 
1) the wave vector making some general angle with the 
symmetry plane of the uniformly distorted hexatic phase. 
The transverse polarization reorientation mode can be 
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described with a simple form fl{p) = cos(kip). The 
direction perpendicular to ki remains liquid. Again for 
definiteness, we present the limit case where ki is taken 
to be ki= ^e^+jBy , a and b being the parameters of 
the "local unit cell" of the phase A. 
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FIG. 2. The 2D section of the phase diagram in coordinates 
d\ and for /i > 0, Eq. ^. Distorted hexatic A, see Fig. 
|l], can give way to crystalline low symmetry phases with con- 
densed transverse polarization modes. In (1) the transverse 
polarization reorientation mode condenses only in one direc- 
tion coinciding with one of the distorted hexatic "unit cell" 
vectors. In phases (2) the angle fl of the polarization reori- 
entation is a sum of two modes with wave vectors out of the 
symmetry plane of the 2D section, where the wave amplitudes 
\fli\ axe in general independent. 

The structure 2b, see Fig. ||, is represented by the 
sum of two transverse polarization modes in the direc- 
tions of ki and k2 with equal amplitudes. The polariza- 
tion reorientation mode thus assumes the form fl{p) = 
fill (cos(kip) -I- cos(k2p)) . Respecting the symmetry 
of the wave vectors we present the limit case where ki 
and k2 are chosen in the form : ki= ^e.j. + jey and 
k2= —^ex+fey . Note that the structures 1, 2a and 
2b can be described using only one angular amplitude 
(see Fig. 2). 

The structure 2c, see Fig. ||, presents the most general 
structure of the ordered phase with two different direc- 
tions of the transverse polarization modes, ki and k2 but 
also with different amplitudes of the waves |r2i| ^ \^2\- 
To present the limit structure we choose the wave vectors 
ki and k2 in the same form as for the previous phase. The 
symmetry is broken with respect to the phases 2a and 
2b : the directions of the transverse polarizations of the 
molecules now make three different angles with respect 
to the y-axis + |n2|, - \^2\ and - + |J72|- 

The phase diagram corresponding to the free energy 



Eq. ^ has the form as presented on Fig. ^ All the transi- 
tions between the distorted hexatic phase and modulated 
phases, as well as between the modulated phases are of 
second order. 

The existence of angular interactions among stiff, ne- 
matically ordered, polymer molecules such as DNA, that 
depend on the orientation of their transverse polarization 
vectors can introduce important modifications into the 
phase diagram at intermediate and high densities. If the 
equilibrium transverse polarization order is homogeneous 
(intermediate densities), corresponding to aligned trans- 
verse polarization vectors on neighboring molecules, than 
the hexatic phase becomes distorted with one axis of the 
local "unit cell" pointing in the direction of the trans- 
verse polarization order. This distortion, that should 
grow as one approaches the hexatic - crystal transition, 
could show up through a broader positional peak in X- 
ray scattering, corresponding to domains of different dis- 
torted hexatic directions, giving a very satisfactory ex- 
planation for observations in DNA 

If on the other hand, the equilibrium transverse polar- 
ization order is inhomogeneous (higher densities) , corre- 
sponding to non-zero angles between neighboring trans- 
verse polarization vectors, the corresponding lattice be- 
comes crystalline with long range positional order either 
in one or both directions and with a deformed hexagonal 
unit cell. Thus instead of having six equivalent nearest 
neighbors, the distorted hexagonal unit cell can give way 
to four nearest neighbors at a non-zero angle (energet- 
ically more favourable) between ils, whereas the more 
distant molecules can have the same direction of ft (en- 
ergetically less favourable). This case too could be as- 
sociated directly with observations in crystallised DNA 
arrays Q . We are thus able to provide a consistent inter- 
pretation for DNA phase behavior in a range of densities 
between the crystalline and the cholesteric phases. 
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